Planck 2018 results. IV. Diffuse component separation

We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow closely those described in earlier papers, adopting four methods (Commander, NILC, SEVEM, and SMICA) to extract the CMB component, as well as three methods (Commander, GNILC, and SMICA) to extract astrophysical components. Our revised CMB temperature maps agree with corresponding products in the Planck 2015 delivery, whereas the polarization maps exhibit significantly lower large-scale power, reflecting the improved data processing described in companion papers; however, the noise properties of the resulting data products are complicated, and the best available end-to-end simulations exhibit relative biases with respect to the data at the few percent level. Using these maps, we are for the first time able to fit the spectral index of thermal dust independently over 3 degree regions. We derive a conservative estimate of the mean spectral index of polarized thermal dust emission of beta_d = 1.55 +/- 0.05, where the uncertainty marginalizes both over all known systematic uncertainties and different estimation techniques. For polarized synchrotron emission, we find a mean spectral index of beta_s = -3.1 +/- 0.1, consistent with previously reported measurements. We note that the current data processing does not allow for construction of unbiased single-bolometer maps, and this limits our ability to extract CO emission and correlated components. The foreground results for intensity derived in this paper therefore do not supersede corresponding Planck 2015 products. For polarization the new results supersede the corresponding 2015 products in all respects

Planck 2018 results. III. High Frequency Instrument data processing and frequency maps

This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous 2015 release. They enabled the first significant measurement of the reionization optical depth parameter using HFI data. This paper presents an extensive analysis of systematic effects, including the use of simulations to facilitate their removal and characterize the residuals. The polarized data, which presented a number of known problems in the 2015 Planck release, are very significantly improved. Calibration, based on the CMB dipole, is now extremely accurate and in the frequency range 100 to 353 GHz reduces intensity-to-polarization leakage caused by calibration mismatch. The Solar dipole direction has been determined in the three lowest HFI frequency channels to within one arc minute, and its amplitude has an absolute uncertainty smaller than $0.35\mu$K, an accuracy of order $10^{-4}$. This is a major legacy from the HFI for future CMB experiments. The removal of bandpass leakage has been improved by extracting the bandpass-mismatch coefficients for each detector as part of the mapmaking process; these values in turn improve the intensity maps. This is a major change in the philosophy of "frequency maps", which are now computed from single detector data, all adjusted to the same average bandpass response for the main foregrounds. Simulations reproduce very well the relative gain calibration of detectors, as well as drifts within a frequency induced by the residuals of the main systematic effect. Using these simulations, we measure and correct the small frequency calibration bias induced by this systematic effect at the $10^{-4}$ level. There is no detectable sign of a residual calibration bias between the first and second acoustic peaks in the CMB channels, at the $10^{-3}$ level

Planck 2018 results. XII. Galactic astrophysics using polarized dust emission

We present 353 GHz full-sky maps of the polarization fraction p, angle \u3c8, and dispersion of angles S of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of p decrease with increasing NH. The uncertainty on the maximum polarization fraction, pmax=22.0% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between p and S is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of p, \u3c8, and S mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map S 7p, looking for residual trends. While p decreases by a factor of 3--4 between NH=1020 cm 122 and NH=2 71022 cm 122, S 7p decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of S 7p with the dust temperature, even though in the diffuse ISM lines of sight with high p and low S tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on NH and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the pmax observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas

Planck 2018 results. XII. Galactic astrophysics using polarized dust emission

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $\psi$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $\psi$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas

Geochemical and Sr-O isotopic constraints on magmatic differentiation at Gede volcanic complex, west Java, Indonesia

The Gede Volcanic Complex (GVC) of the Sunda island arc (West Java, Indonesia) consists of multiple volcanic centres and eruptive groups with complex magmatic histories. We present new petrological, mineralogical, whole-rock major and trace element and Sr–O isotopic data to provide constraints on the relative importance of fractional crystallisation and magma mixing in petrogenesis, as well as on the role and nature of the arc crust. Banded juvenile scoria from Young and Old Gede provide unequivocal evidence for the (late-stage) interaction of distinct magmas at Gede volcano. However, the relatively small-degree compositional zoning observed in plagioclase phenocrysts of all eruptive groups (up to ~20 mol% An) may be attributed to physical changes in magma properties (e.g. P, T, and PH2O) rather than changes in melt composition. Major element and trace element variations within each eruptive series are inconsistent with magmatic evolution through simple mixing processes. Instead, mixing of variably fractionated magma batches is suggested to account for the significant scatter in some element variation diagrams. No correlation is observed between textural complexity and/or mineral disequilibrium and whole-rock geochemistry. REE data and geochemical modelling indicate that fractional crystallisation involving amphibole in the mid- to lower crust, and fractionation of plagioclase, clinopyroxene, Fe–Ti oxide ± olivine ± orthopyroxene provide strong control on the geochemical evolution of GVC rocks. Two-pyroxene geothermobarometry provides pre-eruption crystallisation temperatures of 891–1,046°C and pressures of 3.4–6.5 kbar, equivalent to ~13–24 km depth beneath the volcanoes (mid- to lower crust). Low, mantle-like clinopyroxene δ18O values of GVC lavas and poor correlation of Sr isotope ratios with indices of differentiation precludes significant assimilation of isotopically distinct crust during magmatic differentiation. Therefore, we suggest that the geochemical character of the moderately thick West Javan arc crust is relatively immature compared to typical continental crust. Trace element ratios and strontium isotopes show that the magmatic source composition of the older geographical units, Gegerbentang and Older Quaternary, is distinct from the other GVC groups

Untangling differentiation in arc lavas: constraints from unusual minor and trace element variations at Salak Volcano, Indonesia

Volcanic rocks from Salak Volcano in West Java display intriguing minor and trace element geochemical variations with silica. TiO2 and P2O5 contents, Y, HFSE and REE concentrations are abnormally rich in the Central Vent Group (CVG) lavas (e.g. Y = 32–69 ppm; Yb, 3–6.5 ppm) and display striking positive correlations with SiO2. This contrasts with rocks erupted at side vents (SVG) on the eastern and western flanks of Salak and with rock suites of most other Javan volcanoes where these elements remain relatively constant with increasing SiO2. Modelling of major and trace element data indicate that low pressure fractional crystallisation exerts strong control on the composition of CVG lavas. HFSE and HREE data are inconsistent with magma mixing, and can be explained by incompatible behaviour during fractionation of plagioclase, clinopyroxene and Fe-Ti oxide ± orthopyroxene and olivine. The observed variations in K/Rb and Ba/Th ratios and correlation of 87Sr/86Sr with indices of differentiation necessitate assimilation of a low K/Rb, low Ba/Th, Sr-rich contaminant with 87Sr/86Sr of ~ 0.7048 during fractional crystallisation. For the eastern flank SVG, deep fractionation of a phase in which HFSE and HREE are compatible (e.g. amphibole) is implicated. By extension, this is also suggested to occur beneath the majority of Javan volcanoes. Radiogenic isotope ratios (Sr-Hf-Nd) of Salak lavas are similar to other Javan lavas. SVG rocks erupted from the eastern flank vent have significantly more primitive Hf-Nd isotope ratios than other Salak volcanic rocks and may represent the least contaminated (by arc crust or subduction input) Salak lavas

Constraining fluid and sediment contributions to subduction-related magmatism in Indonesia: Ijen Volcanic Complex

Ijen Volcanic Complex (IVC) in East Java, Indonesia is situated on thickened oceanic crust within the Quaternary volcanic front of the Sunda arc. The 20 km wide calderas complex contains around 22 post-caldera eruptive centres, positioned either around the caldera-rim (CR) or along a roughly NE–SW lineament inside the caldera (IC). The CR and IC lavas exhibit separate differentiation histories. Major element and trace element modelling shows that fractionation of olivine, clinopyroxene, Fe–Ti oxide ± plagioclase is important in the CR group, whereas plagioclase is the dominant fractionating mineral in the same assemblage for the IC group. Spatial controls on magmatic differentiation highlight important structural controls on magma storage and ascent at IVC. Mantle-like 18O values, restricted ranges in Sr, Nd and Hf isotope ratios, and a lack of correlation between isotope ratios and indices of differentiation in IVC lavas are consistent with magmatic evolution through fractional crystallization. Furthermore, the small ranges in isotopic ratios throughout the complex indicate little heterogeneity in the mantle. IVC lavas possess higher 176Hf/177Hf and 143Nd/144Nd isotope ratios than other volcanoes of Java, representing the least contaminated source so far analysed and, therefore, give the best estimate yet of the pre-subduction mantle wedge isotopic composition beneath Java. Trace element and radiogenic isotope data are consistent with a two-stage, three-component petrogenetic model for IVC, whereby an Indian-type mid-ocean ridge basalt (I-MORB)-like fertile mantle wedge is first infiltrated by a small percentage of fluid, sourced from the altered oceanic crust, prior to addition of <1% Indian Ocean sediment dominated by pelagic material